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ABSTRACT 

Current theory of disk galaxy formation is used to study fundamental-plane 
(FP) type of relations for disk galaxies. We examine how the changes in model 
parameters affect these relations and explore the possibility of using such re- 
lations to constrain theoretical models. The distribution of galaxy disks in 
the space of their fundamental properties are predicted to be concentrated in 
a plane, with the Tully-Fisher (TF) relation (a relation between luminosity 
L and maximum rotation velocity V m ) being an almost edge-on view. Us- 
ing rotation velocities at larger radii generally leads to larger TF scatter. In 
searching for a third parameter, we find that both the disk scale-length Rd 
(or surface brightness) and the rotation-curve shape are correlated with the 
TF scatter. The FP relation in the (LogL, LogV^, Log i?d)-space obtained 
from the theory is L ex R°[ Vj^ , with a' ~ 0.50 and j3' ~ 2.60, consistent 
with the preliminary result we obtain from observational data. Among the 
model parameters we probe, variation in any of them can generate signifi- 
cant scatter in the TF relation, but the effects of the spin parameter and 
halo concentration can be reduced significantly by introducing Rd while the 
scatter caused by varying rrid (the ratio between disk mass and halo mass) is 
most effectively reduced by introducing the parameters which describes the 
rotation-curve shape. The TF and FP relations combined should therefore 
provide useful constraints on models of galaxy formation. 

Key words: galaxies: formation-galaxies: structure-galaxies: spiral-galaxies: 
fundamental parameters-dynamics: Tully-Fisher relation 
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1 INTRODUCTION 

Spiral galaxies are characterized by their flattened disks (with approximately exponential 
profiles) and nearly flat rotation curves. The main observational parameters that describe 
the overall properties of a spiral galaxy are the central surface brightness /iq, the disk scale- 
length Rd and the characteristic rotation velocity (e.g. the maximum rotation velocity V m ). 
The observed disk population covers a large range in these parameters, with /x varying from 
~ 24 (low surface brightness galaxies) to ~ 20magarcsec -2 , Rd from ~ 0.1 to ~ lOkpc, and 
V m from ~ 50 to ~ 300 km s -1 . Despite the diversities, spiral galaxies obey a well-defined 
scaling relation between their total luminosity L and maximum rotation velocity V m . This 
relation (called the Tully-Fisher relation, hereafter TF relation) is usually expressed in the 
form 

L = AV? , (1) 

where f3' is the slope, and A is the zero-point. The observed TF relation is quite tight. 
For example, in the /-band, the scatter in luminosities for a fixed V m is only about 0.38 
magnitude (Giovanelli et al. 1997). Thus, the TF relation can be used as a relative distance 
indicator for spiral galaxies. It is still unclear whether the scatter in the observed TF relation 
is dominated by observational errors or by intrinsic variations. Since intrinsic scatter must 
exist to some degree, a natural question is whether there is a third parameter that correlates 
with the residual scatter of the TF relation. The situation is quite similar to that for elliptical 
galaxies. Although elliptical galaxies are found to obey the Faber- Jackson relation (a relation 
between luminosity L and central velocity dispersion a), the fundamental-plane (hereafter 
FP) relation (with the introduction of a third parameter, the effective radius r e ) has much 
smaller scatter (Djorgovski & Davis 1987; Dressier et al. 1987). For disk galaxies, due to the 
broad distribution of scale-length R d , one obvious choice of the third parameter is Rd or the 
central surface brightness /z (/^o an d Rd are related by L = 2n fi R^ for an exponential disk). 
On the other hand, since disk galaxies also present diversities in their rotation-curve shapes 
(Persic, Salucci & Stel 1996), another possible choice may be a parameter that characterizes 
the rotation-curve profile. 

Much effort has been made in searching for this third parameter (e.g. Kodaira 1989; 
Han 1991; Tully & Verheijen 1997; Willick et al. 1997; Courteau & Rix 1999; Willick 1999), 
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but the results are still inconclusive. While most of the studies cited above concluded that 
the third parameter is not crucial, some of them did find evidence for the existence of a 
third parameter (e.g. Kodaira 1989; Willick 1999). More recently, Koda, Sofue & Wada 
(2000a) obtained a scaling plane in the form of L oc (VR) 13 ' from both observational data 
and numerical simulations. But the radius R in this relation is not really a new degree of 
freedom, because it is bounded to the rotation velocity V with the same power index f3'. 

From theoretical side, there are many attempts to understand the formation of disk 
galaxies in the framework provided by current cosmogonic models (e.g. CDM models: Fall 
& Efstathiou 1980; White & Frenk 1991; Dalcanton, Spergel & Summes 1997; Mo, Mao 
& White 1998, hereafter MMW). The current standard scenario of disk formation assumes 
that galaxy disks form as a result of gas cooling in dark matter haloes. Detailed modeling 
shows that such models are quite successful in interpreting observational data, especially the 
observed TF relation (MMW; van den Bosch 1998, 2000; Avila-Reese, Firmani & Hernandez 
1998; Weil, Eke & Efstathiou 1998; Heavens & Jimenez 1999; Mo & Mao 2000; Navarro & 
Steinmetz 2000; Koda, Sofue & Wada 2000b; Zou & Han 2000; Cole et al. 2000; Buchalter, 
Jimenez & Kamionkowski 2001). In the CDM cosmogonies, the characteristic rotation ve- 
locity of a disk galaxy is determined largely by the potential of the dark matter while the 
luminosity of the disk is given by the amount of stars (baryons) in the disk. The observed TF 
relation therefore implies a tight relation between the disk mass and the depth of the dark 
halo potential well, and its scatter may be used to constrain models of galaxy formation. 

Our main goal in this paper is to examine whether FP-like relations exist for the disk 
population based on current theoretical models of disk formation, and to provide some 
theoretical guidelines for the search of FP-like relations. Moreover, we explore the possibility 
of using such relations as a tool to understand galaxy formation by examining the response 
of the theoretical scaling relations to the changes of model parameters. 

The structure of the paper is as follows. Section 2 presents the model and summarizes 
the model parameters. In Section 3, we search for the theoretical FP relation and examine 
its responses to the variations of model parameters. In Section 4, some preliminary results 
obtained from observational data are presented and compared with our model predictions. 
Finally, in Section 5, we make further discussions and summarize our main conclusions. 

2 THEORETICAL MODEL 
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2.1 Disk Formation 

Our model here follows that presented in Mo, Mao & White (1998, hereafter MMW). We 
refer the reader to that paper for details; here we only introduce the main ingredients relevant 
to our analyses. 

The initial density profiles of dark haloes are assumed to take the universal form, 

o( r ) = P critS ° ( 2 ) 

' {) (r/r s )(l + r/r s ) 2 ' {) 

where r s is a characteristic radius, p cr u is the critical density and Sq is a constant (Navarro, 

Frenk & White 1997, hereafter NFW). We define the radius of a dark halo to be r2oo within 

which the mean density is 200 times the critical density, and represent the characteristic 

radius by the concentration, c = r 2 oo/r s . It is then easy to show that 

200 c 3 

3 ln(l + c) - c/(l + c) ' 1 ' 

The radius r 2 oo and mass M of a halo can be expressed in terms of its circular velocity V c as 

V c ,, V c 2 r 2 oo , s 

r20 ° = 10i^) ' M= ~G~ ' (4) 

where G is the gravitational constant, and H(z) is the Hubble constant at redshift z. We 
relate the initial angular momentum J of a halo to its spin parameter A through the definition 

A = J\E\ l ' 2 G- l M-">l 2 , (5) 

where E is the total energy of the halo. A dark halo is therefore described by three param- 
eters: the circular velocity V c , the concentration c and the spin parameter A. As a result 
of dissipative and radiative processes, the gas component gradually settles into a disk. We 
assume the disk mass to be a fraction of the halo mass, and the disk angular momentum 
to be jd times J. Thus, two other parameters, jd and rrid are introduced to relate the halo 
properties to the properties of the final disk. Following MMW, we also assume that disks 
have exponential surface density profiles with a constant mass-to-light ratio T, and that 
dark haloes respond to disk growth adiabatically. Under these assumptions, we can obtain 
the luminosity L, the disk scale length R d and the rotation curve V(R) for a given set of 
model parameters. Specifically, 
L = mdM ^ = 1 / jd_\ = 
T V2 \m d J 

where f r and f v are factors which depend on halo profile and disk action. We use the same 
procedure outlined in MMW to calculate f r and f v as functions of c, jd, m d and A. 
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2.2 Model Parameters 

As discussed above, five parameters (V c , A, c, and jd) need to be specified for individual 
haloes in order to predict the properties of the disks that form within them. In this subsection 
we describe one by one how these parameters are chosen in our modelling. 

2.2.1 Halo circular velocity V c 

The distribution of halo masses is derived from the Press-Schechter formalism (Press & 
Schechter 1974, hereafter PS) applied to the ACDM cosmogony (with mass density param- 
eter Q Q = 0.3, cosmological constant Q\ = 0.7, Hubble's constant h = 0.7, perturbation- 
spectrum normalization erg = 1). We select the host haloes of disks with circular velocities in 
the range from 50kms _1 to 300 km s" 1 at redshift z = 0. The choice of cosmogony does not 
change our main results, because we are not interested in the exact distribution of galaxies 
with respect to mass. 



2.2.2 Halo spin parameter X 

N-body simulations show that the distribution of halo spin parameter A can be approximated 
by a log-normal function, 

/ w j\ 1 [ ln 2 (A/A)l dX 
p(X)dX=^= exp ^J-L — , (7) 

with A=0.05 and o~\ n \ = 0.5. This distribution is found to be quite independent of cosmology 
and of halo mass (Warren et al. 1992; Lemson & Kauffmann 1999). Syer, Mao & Mo (1999) 
obtained a similar distribution function for A from the observational data of disk sizes. In 
this paper, we use equation (^) for the distribution of A. 



2.2.3 Halo concentration c 

The halo concentration, c, may depend on halo formation history and cosmogony. Sim- 
ulations show that the circular-velocity dependence of the NFW concentration c can be 
approximated by 

cnfw ~ cnfw* (K/100 km s" 1 ) ~ 1/3 (8) 

(see NFW). Using iV-body simulations, Jing (2000) found that the distribution of c for the 
majority of his simulated haloes can be described by a log-normal function with a mean c 
slightly smaller than that of the NFW result: 
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p(c')dc' 



1 



In 2 (c'/c' 



dine , c = cj cnfw , 



(9) 




where d = 0.85 and o~\ nc > = 0.25. In this paper, we take cnfw* = 10 as is given in the original 
NFW paper, use equation (||) to obtain cnfw and equation @ for the distribution of c. 

2.2.4 The mass ratio rrid 

This parameter represents the ratio between the disk mass and the total halo mass. The 
processes by which gas cools into a disk to form stars are complicated and not well understood 
at the present time. Nevertheless, some limits can be set. First, the value of rrid should be 
smaller than the overall baryon fraction in the universe f2^/f2 , which is about 0.1 for a 
low-density universe with Qq = 0.3 and h = 0.7, and about 0.05 for an Einstein-de Sitter 
universe with h = 0.5, according to the cosmic nucleosynthesis theory. Second, the value of 
rrid should not be much smaller than 0.01, because disks with lower gas content may not 
be able to form sufficient amount of stars to be included in a TF sample. Thus, we set the 
range of rrid to be from 0.01 to 0.1 and take 0.05 as the typical value. We will also consider 
a case where rrid changes systematically with V c . 

2.2.5 The angular-momentum ratio jd 

This parameter characterizes the ratio between the specific angular momentum of the disk 
and that of the dark halo. Although disk material may have the same initial specific angular 
momentum as the dark matter, little is known about the transfer of angular momentum be- 
tween gas and dark matter in subsequent evolutions. From earlier theoretical considerations 
(e.g. MMW) we know that jd ~ rrid is required in order to ensure the predicted disk sizes 
to match observations. In our model, we allow jd to have 50% chance to be in the range 
0.5md-lrrid, and another 50% chance to be in the range lmd-2rrid. Thus, we adopt for jd a 
simple distribution function of the form: 



2.3 Monte-Carlo Realization 

In order to search for a FP relation for theoretical disks, we need to choose three observ- 
able parameters to define the plane. We take the luminosity L and a characteristic rotation 




djd/md for 0.5 m d < j d < m d 
djd/{2m d ) for m d < jd < 2 m d 



(10) 
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velocity V as the two fundamental quantities. In most of our discussions, we use the max- 
imum velocity V m as the characteristic velocity, but we also discuss the effects of using 
rotation velocities at other radii. For the third quantity, we try two different choices based 
on the discussion in the previous section: the disk scale-length R d and a rotation-curve shape 
parameter V (to be defined below). 

The procedure for obtaining the FP is as follows. We first generate a Monte-Carlo sample 
of 1000 haloes from the PS formalism, with V c in the range between 50 and 300 km s" 1 . We 
then assign to each halo the other four parameters by Monte-Carlo simulations. Finally we 
use the model of MMW to predict the properties of each model galaxy. As in MMW, we 
exclude unstable disks with 

e m = Vm 1/9 <0.9, (11) 

(GM d /R d ) 1/2 

where M d = m d M is the gas mass of the disk (e.g. Christodoulou, Shlosman & Tohline 
1995). On the other hand, galaxies with very low density may not form stars. We therefore 
also exclude disks with Toomre parameter 

«^>1.4, (12, 

where a is the gas velocity dispersion, k is the epicyclic frequency, and S is the mass surface 
density (Kennicutt 1989). In our model, we use the average surface mass density inside the 
half-mass radius (about 1.68R d ) to replace S, and adopt a = 6kms _1 . We assume that all 
disks have the same mass-to-light ratio and take T l.7h in the /-band (Bottema 1997). 
All the results are for disks at z = and we take h = 0.7 whenever an explicit value of h is 
used. 



3 MODEL PREDICTIONS 

3.1 Results Based on Maximum Rotation Velocity and Disk Scale-length 

In this section, we try to find the FP of disk galaxies in the ( Log L- Log V m - Log R d )-space. 
The reason is clear: while L and V m are the two parameters in the conventional TF relation, 
the choice of R d as the third parameter is motivated by the fact that disk galaxies cover a 
large range of R d . We define the FP as 

Mj = a Log R d ( kpc) + Log V m ( km s" 1 ) + 7, (13) 
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Table 1. The effects of changing individual model parameters on the TF relation and the FP relation defined by equation 
jl3). Column 1 lists the name of each case, columns 2 to 6 give the model parameters (which are either a single number, or a 
range, or an equation of distribution), columns 7 to 9 show the fitting results for the FP relation, while columns 10 and 11 list 
the scatter of the magnitude Mj in the FP and TF relations, respectively 



Case 


A 


3d 


c 


m d 


T 


a 


P 


7 


crpp 


(TTP 


A 


eq. (7) 


0.05 


8.5 


0.05 


1.7/i 


-1.29 


-6.16 


-7.21 


0.08 


0.39 


3d 


0.05 


eq. (10) 


0.05 


8.5 


1.7/i 


-1.21 


-6.27 


-7.05 


0.04 


0.29 


c 


0.05 


0.05 


eq. (9) 


0.05 


1.7/i 


-2.89 


-4.61 


-9.83 


0.002 


0.23 


m d 


0.05 


m d 


8.5 


0.01-0.1 


1.7/i 


1.11 


-8.53 


-3.23 


0.12 


0.20 



where Mi is the absolute magnitude of the disk in the /-band, a, j3 and 7 are constants 



which are obtained by a least-square fit of equation (|13| ) to the Monte-Carlo data. In the 
disk model outlined above, there are four parameters, A, m^, and c, whose effects on 
the FP relation need to be evaluated. If all these four parameters are set to be constant, 
the predicted TF relation is a line with no scatter. The scatter in any one of these four 
parameters can cause scatter in the TF relation (see e.g. MMW, Mo & Mao 2000). The first 
question is whether some of the scatter can be eliminated or reduced in the FP relation. In 
order to show this, we first construct four samples where scatter is allowed only in one of 
the four model parameters while the other parameters taking their typical values (see Table 
1). In these cases, we ignore the selection criteria given by equations (|TT| ) and fll2|), and so 
some of the model galaxies in these samples may not correspond to any realistic disks. The 
scaling relations are shown in Figure 1, while the fitting results of these relations are listed 
in Table 1. The scatters on both the TF relation and the FP relation are represented by the 
root-mean-square deviation (rms) of Mj. 

Comparing the scatter of the TF relation and the corresponding FP relation, we see that 
the introduction of Rj as the third parameter is very effective in reducing the scatter caused 
by A, jd and c but less so in reducing the scatter caused by rrid- Notice that the FPs in 
different cases are tilted significantly with each other, except cases of varying A and jd (see 
Table 1). 

After an preview of the effects of individual model parameters on the scaling relations, 
we now examine cases where scatter is allowed in more than one model parameters. For 
clarity, we list all cases in Table 2. In all these cases, the selection criteria given by equations 



PI ) and (|T2"D are imposed. 
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Figure 1. The effects of changing individual model parameters (as indicated in the panels) on the TF and FP relations, (a), 
(b), (c) and (d) show the effects of A, j d , c and m d respectively. In each case, the left panel shows the FP relation (the magnitude 
Mi as a function of aLogRj + /3Log V m + 7, where the values of a, (3 and 7 arc listed in Table 1), whereas the right panel 
shows the corresponding TF relation (the solid line being the linear regression). The scatter in Mi is represented by the rms 
deviation. 



Table 2. The predicted TF (Mr-Log V m ) and FP (Mj- Log V m - Log R d ) relations in various cases (Case VI will be discussed 
in Subsection 3.4). The table structure is the same as Table 1. 



Case 


A 


3d 


c 


m d 


T 


a 


P 


7 


OTP 


CTTF 


I 


cq. (7) 


eq. (10) 


8.5 


0.05 


1.7ft 


-0.95 


-6.54 


-6.64 


0.03 


0.21 


II 


cq. (7) 


cq. (10) 


eq. (9) 


0.05 


1.7/i 


-1.22 


-6.43 


-6.65 


0.15 


0.35 


III 


cq. (7) 


cq. (10) 


8.5 


0.01-0.1 


1.7/i 


-1.27 


-6.35 


-6.79 


0.32 


0.44 


IV 


cq. (7) 


cq. (10) 


eq. (9) 


0.01-0.1 


1.7/i 


-1.48 


-6.35 


-6.61 


0.34 


0.50 


V 


cq. (7) 


cq. (10) 


cq. (9) 


0.01-0.1 


cq. (15) 


-0.79 


-6.47 


-6.85 


0.33 


0.38 


VI 


cq. (7) 


eq. (10) 


cq. (9) 


eq. (18) 


1.7/i 


-1.09 


-8.11 


-2.90 


0.22 


0.35 



3.1.1 Effects of X and jd 



Since the TF relation is, for a constant disk mass-to-light ratio, a relation between the disk 
mass and the maximum rotation velocity, difference in the disk angular momenta may cause 
scatter in the TF relation because, for a given disk mass, the contribution by the disk to the 
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Figure 2. The predicted TF and FP relations for Case I in Table 2. The predicted TF relation is shown in panel (a) together 
with the fit to the model prediction (solid line) and the fit to the observational result of Dale et al. (1999b) (dashed line). 
Panel (b) shows the FP with as the third parameter, while panel (c) and (d) show the FPs based on the two different 
rotation-curve shape parameters Ti and Y2 (to be discussed in Subsection 3.3). 

rotation velocity depends on disk angular momentum. The disk specific angular momentum 
is characterized by the product Xja and so we examine the effects of these two parameters 
together. We generate a Monte-Carlo sample, with A and following the distributions 
described in equations (0) and fllPf) , while keeping = 0.05 and c = 8.5 (see Case I in 



Table 2). 

Panel (b) of Figure 2 shows the FP, while the corresponding TF relation is shown in panel 
(a) together with a comparison to the observational result of Dale et al. (1999b). We see that 
there is an almost perfect FP in this case, although there is significant TF scatter caused 
by the variations of the spin parameter A and angular-momentum ratio jd- The observed 
slope and zero-point of the TF relation are well reproduced, which is consistent with the 
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Figure 3. The predicted FP and TF relations for Case II in Table 2. The setting of the panels is the same as in Figure 2. 

result already found in MMW. The FP is given by a = —0.95, (3 = —6.54, with a zero-point 
7 = —6.64. The FP scatter in Mj is only 0.03 mag, much smaller than the predicted TF 
scatter (0.21 mag). Thus, the scatter in the TF relation caused by Xj^ can be eliminated 
almost entirely by the third parameter Rd- So, if the scatter in the TF relation were caused 
entirely by the dispersion in the spin, the model would predict a perfect FP for the disk 
population. However, the predicted TF scatter (0.21 mag ) is lower than the observed value 
of Dale et al. (0.38 mag, see also Giovanelli et al. 1997), implying that other factors are also 
important in causing the observed TF scatter. 

3.1.2 Effect of c 

Here we include the scatter caused by the distribution of c, in addition to those by A and ja, 
but still keep = 0.05. The results are shown in Figure 3. In this case (Case II in Table 
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Figure 4. The predicted FP and TF relations for Case III in Table 2. The setting of the panels is the same as in Figure 2. 

2), the observed TF relation is also well reproduced. The predicted FP (with a = —1.22, 
(3 = —6.43 and 7 = —6.65) is close to Case I shown in Figure 1, except the scatter is larger 
(0.15 mag). This scatter in the FP is still substantially smaller than that in the predicted 
TF relation (0.35 mag). The introduction of Rd cannot eliminate all the scatter, because the 
FP defined by c is tilted with respect to that defined by A and jd. ( see Table 1). 

3.1.3 Effect of m d 

To examine the effect of on the scaling relations, we construct a Monte-Carlo sample with 
a fixed concentration (c = 8.5), with A and jd having their typical distributions, and with 
irtd randomly drawn from 0.01 to 0.1 (Case III in Table 2). As discussed in Subsection 2.2.4, 
this range of md covers the values expected from any reasonable considerations. We keep c 
constant, in order to single out the effect of m^. The results are shown in Figure 4. In this case, 
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Figure 5. The predicted FP and TF relations for Case IV in Table 2. The setting of the panels is the same as in Figure 2. 

the slope and zero-point of the FP (a = —1.27, j3 = —6.35 and 7 = —6.79) are comparable 
to those shown in Figures 2 and 3. The scatter in the predicted TF relation (0.44 mag) is 
much larger than those in Case I and II and even larger than that in the observations of 
Giovanelli et al. (1997). This large predicted scatter is obviously a consequence of the large 
range of used in this case. The scatter in the FP relation (0.32 mag) is also quite large, 
suggesting that the scatter produced by the variation of rrid has a large tilt with respect to 
that caused by A and jd (see Table 1). 

3.1.4 A general case 

As a summary, we consider a case where scatter is allowed in all of the model parameters 
(Case IV in Table 2). The results are shown in Figure 5. Although generous amount of 
scatter is assumed for each of the model parameters, a FP can still be defined for the model 
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galaxies, with a = —1.48, /3 = —6.35 and 7 = —6.61. The predicted TF relation has a 
scatter (0.50 mag) larger than that observed, again mainly due to the large range of md 
assumed. The scatter of the FP is also quite large (0.34 mag). This is not surprising, because 
the scatter is dominated by the variation in m d . 

If we express the FP in the form L oc R^'V? , the slopes a' and (3' change only little for 
different cases, and the FP relation for the model disks can be represented by 

L oc B$V£ , with a' = 0.49 ± 0.10 , (3' = 2.60 ± 0.05 , (14) 

where the errors represent the scatter among different cases. Thus, the luminosity of a disk 
galaxy is mainly determined by V m , and so the TF relation is almost an edge-on view of 
the FP. However, the residual dependence on R d is also significant, due to the fact that the 
range of Rd for the observed disks is quite large. 

3.1.5 A case with varying T 

Although we did not include the disk mass-to-light ratio T in the list of our model parameters 
above, it obviously is another parameter in the model. In order to compare model predictions 
directly with observations, we need to convert the predicted disk mass to a disk luminosity, 
which requires an assumption about the disk mass-to-light ratio. Scatter in the disk mass- 
to-light ratio therefore also gives rise to scatter in both the TF and the FP relation. Clearly, 
the induced scatter is along the Mj axis, with a rms (in terms of the luminosity) exactly 
the same as that in the mass-to-light ratio. 

The exact value of T for a galaxy depends on its star formation history and is in general 
difficult to model accurately. As an example, we consider a case (Case V in Table 2) where 
the disk mass-to-light ratio is assumed to be 

T = 1.7/1(7^/0.05) (15) 

and all other parameters are assumed to be the same as in Case IV. In this case, the effect of 
m d on the magnitude Mj is eliminated, and so the TF scatter caused by rrid is only through 
its effect on the rotation curve. The results are shown in Figure 6. The predicted TF relation 
has scatter smaller than that in Case IV because the effect of m d on Mj is eliminated [see 
Equation (6)]. The scatter of the predicted FP (0.33 mag) is close to that of the predicted 
TF relation (0.38 mag), implying that the scatter in the TF relation caused by rrid through 
disk action cannot be effectively reduced by introducing Rd as the third parameter. 

From the results presented above we can conclude that variations in the mass ratio m d 
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Figure 6. The predicted FP and TF relations for Case V in Table 1. The setting of the panels is the same as in Figure 2. 

and in the mass-to-light ratio may be important sources of scatter in the predicted TF and 
FP relations. If the variations in rrtd and T are small, the scatter in the predicted FP relation 
may be significantly lower than that in the TF relation. On the other hand, if the scatter 
in the observed TF relation is mainly due to the variations in and T, the introduction 
of the third parameter will not reduce the scatter significantly. Thus, by comparing the 
scatter in the TF relation with the scatter in the FP relation, one may hope to find the main 
sources of scatter in the two scaling relations. We will come back to this issue in Section 4. 

3.2 Results Based on Rotation Velocities at Other Radii 

Although the maximum rotation velocity (V m ) is commonly used in defining the TF relation 
(mainly for observational reasons, because the inclination-corrected 21 cm line width of a 
disk galaxy is assumed to be twice its maximum rotation velocity), there is no a priori reason 
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Figure 7. The predicted TF and FP relations based on different rotation velocities. The right panels show the FPs with 
as the third parameter, using respectively V1.5, V2.2 and V5 (all in units of kms -1 ) as the characteristic rotation velocity. The 
left panels show the corresponding TF relations. The solid lines in the TF panels are the fits to the model predictions, while 
the dashed lines represent the relation observed by Dale et al. (1999b). 

against using rotation velocities at other radii. Such definition is possible for disk galaxies 
with measured rotation curves. Since the disk contribution to the rotation curve may change 
with radius (for example, the rotation curve of a galaxy at large radius is expected to be 
dominated by dark matter while the disk contribution may be significant in the inner region), 
analyses of the TF relation (or the FP relation) using rotation velocities at different radii 
may provide more information on the mass components of disk galaxies. 

Observationally, there are some attempts to define the TF relation using rotation veloc- 
ities other than the maximum rotation velocity, in the hope of getting smaller scatter. For 
example, Courteau (1997) found an optimal TF relation using rotation velocities at 2.2i?^, 
while there are also claims that the TF relation is optimized if rotation velocities at smaller 
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Table 3. The same as Table 2 except that other rotation velocities are used instead of V r 



Case 


A 


jd 


c 




T 


a 


P 


7 




°TF 


Vl.5 
V 2 .2 

v 5 


cq. (7) 
cq. (7) 
eq. (7) 


eq. (10) 
cq. (10) 
eq. (10) 


eq. (9) 
eq. (9) 
cq. (9) 


0.01-0.1 
0.01-0.1 
0.01-0.1 


1.7h 
1.7h 
1.7/i 


-1.40 
-1.48 
-1.65 


-6.63 
-6.56 
-6.32 


-6.29 
-6.23 
-6.74 


0.29 
0.30 
0.39 


0.44 
0.47 
0.57 



radii are used (see Willick 1999 and references therein). Interestingly, almost all of these 
analyses showed that the scatter in the TF relation becomes larger when the amplitudes of 
the rotation curves at large radii are used. 

In this section, we use our theoretical model to examine the effect of using different 
rotation velocities on the TF and FP relations. Specifically, we choose three representative 
radii, R = 5Rd, 2.2Ra, l.bRd, to define the characteristic rotation velocities (V5, V2.2 and 
V1.5). For most of our model rotation curves, the maximum rotation velocity is reached at 
R ~ 3Rd (MMW), and so V5 represents the rotation velocity well beyond the peak of a 
rotation curve. V2.2 is chosen because the rotation velocity of a pure exponential disk peaks 
at R m 2.1QR d , while V1.5 is chosen to represent the inner part of the rotation curve. 

We construct a Monte-Carlo sample in the same way as described in Subsection 2.3. The 
distributions of all the model parameters are assumed to be the same as the general case 
(Case IV) in Table 2. The results are shown in Figure 7 and summarized in Table 3. 

As shown in Figure 7, the model predictions are consistent with the observations that 
the use of rotation velocities at larger radii generally leads to larger scatters in the TF 
relations. This result can be understood as follows. At very large radii, the disk contribution 
to the rotation curve is negligible, and so the variation in does not affect the rotation 
velocities very much. Since the luminosity of a disk is directly proportional to (for a 
constant mass-to-light ratio), the variation in m d affects the TF scatter directly. This effect 
is reduced in the cases where rotation velocities in the inner part (e.g. V2.2 and V1.5) are 
used, because of the increasing halo-disk interaction (Mo & Mao 2000; Navarro & Steinmetz 
2000). We have also made calculations for models where m d is kept constant. In such cases, 
the TF scatter is almost independent of the definition of the characteristic rotation velocity. 
Thus, the scatter in the TF relations defined at different rotation velocities may be used to 
constrain the variation of the mass ratio m<j. 

Similar to the TF relations, the scatter in the FP relation is also larger when the rotation 
velocities at larger radii are used. This is because R d cannot effectively reduce the TF scatter 
caused by m^. 
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Table 4. The predicted FP relations where the rotation-curve shape parameters Ti and T2 are adopted as the third parameter. 
Column 1 lists the name of each case as listed in Table 2. Columns 2 to 5 give the fit results for the FP in the Mj- Log V m Ti 
space, while columns 6 to 9 are those for the FP in the Mj- Log space. 



Case 



TP 



crpp 



I 


6.76 


-7.31 


-6.16 


0.19 


-9.00 


-7.38 


3.19 


0.19 


11 


-1.05 


-7.55 


-4.93 


0.35 


-14.87 


-7.05 


9.13 


0.26 


III 


10.57 


-7.33 


-6.44 


0.21 


8.07 


-7.32 


-13.31 


0.27 


IV 


8.84 


-8.05 


-4.68 


0.39 


6.71 


-7.97 


-10.55 


0.41 


V 


-5.56 


-7.62 


-4.41 


0.32 


-6.43 


-7.74 


1.56 


0.26 


VI 


-0.28 


-9.97 


0.44 


0.35 


-8.57 


-10.69 


1.94 


0.34 



3.3 Rotation-Curve Shape as the Third Parameter 

Since rotation curves of individual galaxies have different shapes, it is possible that the TF 
scatter is correlated with the rot at ion- curve shape. Persic & Salucci (1990) considered this 
possibility from observational data. In the theoretical model we are considering here, the 
shapes of rotation curves can be affected by all the four model parameters, and so the TF 
scatter must be correlated to the rotation-curve shape to some degree. In this subsection, 
we carry out a quantitative analysis of such correlation. To do this, we use two parameters 
similar to those defined by Persic, Salucci & Stel (1996), 

r 1 = (^iZ) . r 2 = ^, (16) 

\dhogRj 22Rd V 2 . 2 
to represent the rotation-curve shapes in the inner and outer regions. We investigate the FP 
in the space spanned by Mi, Log (V m ) and T\ (r 2 ) for all the cases considered above. The 
simulated FP is fitted to the relation 

Mj = ar + /3Log\/ m (kms- 1 ) +7, (17) 

where T is either F 1 or T 2 . The predicted FPs are shown in panel (c) and panel (d) in each 
of Figures 2-6, and the fitting results are summarized in Table 4. 

These results show that the introduction of the shape parameter of rotation curves as 
the third parameter can reduce the scatter of the TF relation, but the two parameters 
considered here are less effective than Rd in all the cases except in Case III where the TF 
scatter is effectively reduced by the introduction of Ti or T2- This suggests that these shape 
parameters are effective in reducing the scatter caused by m d (due to disk action) but not 
so much in reducing the scatter caused by Xjd and c. 

Unfortunately, Ti and T 2 are both more difficult to measure than R d from observations, 
and it is not yet possible to obtain such a FP from observational data. Notice that the 
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Figure 8. The edge-on view of the FP in the space spanned by Mj, Log-R^i Log V m , r (or for Case IV. Notice that the 
scatter in Mj is reduced significantly relative to that shown in Figure 5. 

slope a associated with the two shape parameters is very sensitive to the change of model 
parameters, because the ranges of these parameters are relatively small. 

From the above discussion we see that the disk scale-length R d is effective in reducing 
the TF scatter caused by Xjd (Figure 2) and c (Figure 3), while the rotation-curve shape 
parameters are more effective in reducing the scatter caused by nid (Figure 4). It is therefore 
interesting to see what happens if both Rd and one of the rotation-curve shape parameters 
are introduced to define a plane in the four-dimensional space spanned by Mi, Log Rd, 
LogV^, and T 1 (or T 2 ). The results are shown in Figure 8 for Case IV. As expected, the 
scatter is significantly reduced. 
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Figure 9. The predicted FP and TF relations for Case VI in Table 2. The setting of the panels is the same as Figure 2 

3.4 A case where changes with V c 

In reality, the value of may depend systematically on halo circular velocity, as is the case 
if feedback effects can eject gas more effectively from smaller haloes (e.g. White & Frenk 
1991). To model this effect, we construct a sample where is assumed to change with V c 
as 

-l 

(1 + 5), (18) 

where 5 is a random number between —0.1 and 0.1. Other parameters are chosen the same 
as in Case IV. This case is listed as Case VI in Tables 2 and 4, and the results are shown 
in Figure 9. In this case the TF relation is Mj = — 9. 75 Log V m — 0.03 [the solid line in 
panel (a) of Figure 9], which has a much steeper slope (f3 ~ 3.8) than all other cases we have 
considered. Correspondingly, the slope (3 in the FP relation (with Rd as the third parameter) 
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Figure 10. The correlation between Mj and the rotation-curve shape parameters, Ti and I^. The top two panels are for Case 
IV while the lower two panes are for Case VI where is correlated with V c as in Equation (18). 

is also steeper, but the value of a does not change much (see Table 2). The introduction of 
Ti (or r 2 ) as the third parameter does not help much in reducing the TF scatter because, 
with the assumed small scatter in rrid, the TF scatter is dominated by the variations in \jd 
and c. 

Since the rot at ion- curve shape parameters are correlated with rrid, it is interesting to see 
if these parameters can be used to reveal the trend of rrid with V c . To do this, we examine 
the correlation between Mj and T 1 (or T 2 ) for Case VI and compare the results with those 
for Case IV. The results are shown in Figure 10. As one can see, if m d is correlated with 
V c , then there is a strong correlation between Mj and Ti (or r 2 ), in the sense that more 
luminous galaxies have smaller T\ (or r 2 ). Such correlations are in fact observed by Persic, 
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Figure 11. The observed and simulated FP and TF relations based on the data from Dale et al. (1999a). Panel (a) shows the 
FP while panel (b) shows the TF relation. Panels (c), (d), (e) and (f) show the FPs and TF relations for two simulated samples 
which have the same distributions in L and Rd as the observed sample. Panels (c) and (d) show the theoretical predictions of 
Case IV. Panels (e) and (f) show the theoretical predictions of a sample with the same parameter as Case II except that the 
mass- to-light ratio T is assumed to be random between Ih and 2h. The solid lines in three TF panels are the direct fit to the 
data point, while the dashed lines show the result obtained by Dale et al. (1999b) after bias correction. 

Salucci & Stel (1996), but the data are still too uncertain to give any meaningful constraints 
on the model. 



4 COMPARISON WITH PRELIMINARY OBSERVATIONAL RESULTS 

Having seen the theoretical expectations for the scaling relations of spiral galaxies, one is 
obviously interested in seeing what observations actually show. Here we present some pre- 
liminary results without going into the details of the observational selection effects, leaving 
a more sophisticated analysis to a future paper. 

The data we use are from Dale et al. (1999a), which presents TF observations for 35 rich 
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Abell clusters of galaxies in the /-band, along with the exponential scale-lengths Rd- We 
simply use the published data without correcting for any bias. Panel (a) and (b) of Figure 
|Tl~| show the results. The FP for the observed galaxies can be described as 

Mi = -1.01 Log i? d (kpc) -5.89 Log W( km s" 1 ) -6.49, (19) 

where W is the velocity width, which we assume to be twice V m . The index on Rd is 
comparable to that given by the models presented above, but the index on W (or V m ) is 
lower than the model predictions. This lower value is almost certainly due to observational 
bias. In order to show this, we construct a sample with the same model parameters as in 
the general case (Case IV) and select model galaxies with the same luminosity and size 
distributions as the observed sample so as to reproduce the selection effect. The results are 
shown in panels (c) and (d) of Figure 11. As one can see, the TF and FP slopes for this 
sample are quite close to the observed values. Note that the scatter is smaller than that 
of Case IV in Figure 5, indicating the selection effect is important. As another example, 
we show in panels (e) and (f) the results for a sample with the same model parameters as 
in Case II except that T has a random distribution between lh and 2h and with galaxies 
also selected according to the observed luminosity and size distributions. We see that the 
scatter allowed in T (a factor of two) is much smaller than that in (a factor of ten), again 
because disk action reduces TF scatter (Mo & Mao 2000). In both cases, the scatter in the 
FP is significantly smaller than that observed. Unfortunately, it is unclear how seriously 
these discrepancies should be taken. The bias correction required may be more complicated 
than what is assumed here. In particular, the observed sample is for cluster galaxies, and so 
systematic bias may also arise from some environmental effects which are not modelled in 
the theory. 

5 DISCUSSION AND SUMMARY 

In this paper, we use current theory of disk galaxy formation to study whether a FP- 
like relation is expected for spiral galaxies. After examining in detail how the changes in 
model parameters affect the scaling relations of disk galaxies, we find that the fundamental 
properties of disk galaxies are generally concentrated into a plane, with the TF relation 
representing an almost edge-on view. We made a systematical search for the third parameter 
which may correlate with the TF scatter. We find that the disk scale-length Rd (or surface 
brightness) as the third parameter can effectively reduce the scatter of TF relation, especially 
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that caused by the variations of halo spin and concentration. The rot at ion- curve shape 
as the third parameter can be used to reduce the scatter caused by rud- For the various 
cases we analyzed, the FPs in the (LogL, LogV m , Logi?d)-space are quite similar and well 
represented by equation (|I~4"D. This relation is consistent with the preliminary results we 
obtain from observational data. 

We also discuss the possibility of using other characteristic velocities [e.g. V5, V2.2 and 
V1.5] instead of V m as the velocity parameter in the TF and FP relations. The theory is 
consistent with the observations that the scatter in the TF relation is generally larger when 
rotation velocities at larger radii are used. 

There are, however, a number of uncertainties in the theoretical model, which must 
be taken into account when comparing models with observations. Real galaxies may be 
much more complicated than our simple model implies: galaxy disks may not be perfectly 
exponential, the mass distribution in dark haloes may be non-spherical and clumpy and so 
the rotation curves may not be smooth, and the existence of galactic bulges may affect both 
the rotation curve and the luminosity profile. Furthermore, the assumption of a constant 
disk mass-to- light ratio (in the /-band), although consistent with current observations, is 
clearly unrealistic, because the mass-to-light ratio of a galaxy depends on its star formation 
history. None of these issues are easy to resolve, but we hope the present paper can provide 
some theoretical guidelines for the search of scaling relations for disk galaxies. 
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